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Abstract. The temporal evolution of a water-sand interface driven by gravity is experimentally investi¬ 
gated. By means of a Fourier analysis of the evolving interface the growth rates are determined for the 
different modes appearing in the developing front. To model the observed behavior we apply the idea of 
the Rayleigh-Taylor instability for two stratified fluids. Carrying out a linear stability analysis we calculate 
the growth rates from the corresponding dispersion relations for finite and infinite cell sizes. Based on the 
theoretical results the viscosity of the suspension is estimated to be approximately 100 times higher than 
that of pure water, in agreement with other experimental findings. 

PACS. 47.54.-l-r Pattern selection - 47.55.Kf Multiphase and particle-laden flows - 68.10.-m Fluid surfaces 
and fluid-fluid interfaces - 81.05.Rm Granular materials - 83.70.Hq Suspensions 


1 Introduction 

There has been great interest in the behavior of granu¬ 
lar materials over the last years (for a review see Q and 
references therein). Examples for the surprising behavior 
of granular matter are two-dimensional localized states, 
called oscillons |^, in vertically vibrated containers, strat¬ 
ification phenomena observed while pouring granular mix¬ 
tures onto a pile or singing sand audible in deserts Q . 

In contrast to the amount of phenomena one lacks a 
sound theoretical explanation for the observations. Diffi¬ 
culties arise due to the highly complex, disorded structure 
formed by the grains and their nonlinear internal friction. 
Therefore, most theoretical approaches have been done by 
numerical methods. Molecular-dynamics simulations and 
cellular automata calculations are frequently applied, for 
recent reviews see . These methods describe in some de¬ 
tail the interactions between the particular grains . Con¬ 
clusions on the macroscopic behavior of granular matter 
are then drawn from the simulations with a great num¬ 
ber of particles. Despite the recent progress in computer 
performance the size of the systems is still too small for a 
quantitative comparison with real experiments 

Hydrodynamic approaches to granular material are few 
and are associated with restrictions as no in¬ 
terparticle correlations B , a Gaussian distribution in the 
velocity of the grains or steady-state properties |pl| . 
Nevertheless, there are striking phenomenological similar¬ 
ities in the observed patterns for pure granular materi¬ 
als and pure fluids. Experiments with an inclined chute 
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[^[l^ or a vertically vibrated container ||^,|^ show the 
most noteable analogy. It is clear that granular media are 
different from fluids but under certain conditions these 
differences are not prevailing. 

The flow of grains embedded in an interstitial fluid 
is dominated either by the effects of grain inertia or by 
effects of fluid viscosity. The Bagnold number B [|T6| ex¬ 
presses the ratio of collision forces between the grains to 
viscous forces in the fluid-grain mixture. A small Bagnold 
number, B < 40, characterizes the regime of the macro- 
viscous flow. In this regime the viscous interaction with 
the pure fluid is important. Examples for this type of flow 
are mud slides and the transport of water-sand mixtures 
in river beds. At large Bagnold numbers, B > 450, the 
flow is called grain-inertia regime where the grain-grain 
interactions dominate. All flows of grains with air as in¬ 
terstitial fluid fall into the grain-grain regime. 

Here, experiments are performed with sand dispersed 
in water. The occuring shear rates, the mean particle di¬ 
ameters, and the viscosity of water result in a Bagnold 
number of about 1 |0. This motivates the idea to con¬ 
sider the water-sand mixture as fluid-like. In the experi¬ 
ments we observe that the initial flat water-sand interface 
evolves into a finger-like pattern. The measured velocity 
for the largest finger was typically three times the Stokes 
velocity of a settling hard sphere. Furthermore, the veloc¬ 
ity was nearly independent of the mass of the sand and 
the diameter of the particles ||^. This independence of 
particle properties encourages a fluid-like description of 
the water-sand mixture. To model this behavior we test 
a continuum approach which is based on a well-known 
hydrodynamic instability. 
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The aim is to determine the critical parameters for the 
stability of the pattern and the dispersion relation for the 
whole spectrum of wave numbers detectable in the experi¬ 
ment. Particularly, the latter extends significantly the ob¬ 
jectives of earlier experiments with water and polysterene 
spheres [Q and with water and glass beads Q. The de¬ 
pendence of the initial wavelength of the developing; pat¬ 
tern on the width of the cell was studied in jl^. In 
the experiments were focused on the evolution of voidage 
shock fronts caused by a step increase or decrease in the 
fluidization velocity for a fluidized bed. We, however, turn 
our attention to the temporal evolution of all wave num¬ 
bers. 

To model the water-sand mixture as a Newtonian fluid 
with effective properties depending on the concentration 
of the particles is obviously a simplifying description. In 
such a model it is assumed that the particle concentra¬ 
tion in the flowing mixture is almost constant. Further¬ 
more, the particles have to be large enough to neglect their 
Brownian diffusion. Stability analyses with these simpli- 
hcations were made for moderate concentrated mixtures 
in horizontal Hagen-Poiseuille flow ||^ and in sedimenta¬ 
tion problems in inclined narrow channels [^l],|2|,^. For 
the latter configuration even the sediment was described 
by such a simple model and a fair agreement with exper¬ 
imental results was found |Q. In the clear knowledge of 
the limitations of a fluid model for the water-sand mixture 
we will examine whether such an approach can catch the 
essence of the experimental results. 

In the following section, the experimental arrangement 
is described and the results for the growth rates of the 
wave numbers are presented. In Sect. ^ the model is ex¬ 
plained and thereafter the calculated growth rates are 
compared with those of the experiment (Sect. The final 
section contains our conclusions and some remarks about 
further prospectives. 

2 Experiment 

2.1 Experimental setup 

A closed Hele-Shaw-like cell is used to investigate the tem¬ 
poral evolution of a water-sand interface driven by gravity 
(see Fig. |^. The cell, a CCD-camera and a neon tube are 
fixed to a frame which can be turned around a horizontal 
axis. This allows image analysis in the comoving frame 
and ensures a homogeneous illumination at every stage of 
the pattern forming process. The length of the cell is 160 
mm, the height 80 mm and the width 4 mm, respectively. 
The cell is filled with sand and destilled water. As sand 
we use spherical glass particles (Wiirth Ballotini MGL) 
of different sizes and size distributions (see Table for 
details). Its material density is given by 2.45 g/cm^. The 
rotation axes of the frame is right beneath the sand layer. 
This minimizes the centrifugal forces on the sand layer 
while the cell is turned. The cell is rotated by hand. To 
obtain reproducible results the vertical and horizontal ac¬ 
celeration is measured by acceleration sensors (ADXL05). 
When the cell passes a rotation angle of 170 degree, a 


lightgate triggers a number of snapshots. This moment 
defines the starting time of our measurements where we 
take images every 20 ms for later analysis. The images 
have a dimension of 256 x 300 pixel. To achieve a reason¬ 
able resolution we only focus on a horizontal length of 61 
mm at the middle of the cell. This gives a resolution of 4.9 
pixel/mm. After each rotation the suspension sediments 
until all particles are settled. This process takes typically 
less than 1 minute. By means of tracer particles it was 
proven that at this time the fluid in the cell is at rest. 
The time between consecutive runs is at least 3 minutes 
to secure independent runs. The influence of the waiting 
time on the occuring patterns was tested and no effect 
concerning the wavelength of the fingers was observed. 

2.2 Experimental results 

Fig. H shows typical images of the sand-water cell at cer¬ 
tain stages. 20 ms after the series of snapshots is started, 
the initial flat sand layer is modulated at small scales (Fig. 
1(a)). These disturbances are enhanced and give rise to 
sand fingers as seen in Fig. H(c). At later stages the fin¬ 
gers evolve to a mushroom-like pattern (Fig. ||(d) and Fig. 
y(e)). This type of pattern was found also by numerical 
simulations for two stratihed suspensions of different con¬ 
centration Q. 

To analyze this behavior we apply a threshold algo¬ 
rithm to obtain the water-sand interface. We look at ev¬ 
ery column of our digitized image to determine the point 
where the grey scale exceeds a certain value. We start at 
the bottom (water) and continue to the top (sand). In 
this way we track down the interface of the pattern. Fig. 
^ shows the temporal evolution of the images presented 
in Fig. ||. Here the interfaces of all patterns are shown. 
While our detecting method works for the patterns dur¬ 
ing the first stages, namely small scale modulations and 
sand fingers, it breaks down for mushroom-like patterns. 
However, this is not crucial because those patterns are 
beyond the scope of the linear analysis presented here. 
Discrete-Fourier-Transformation (DFT) gives the Fourier 
spectrum of each interface. Fig. ^ shows the temporal evo¬ 
lution of the amplitude A of a typical Fourier mode. It is 
seen that A grows exponentially from the first image to t 
= 200 ms. By an exponential fit 

A{k) = Ai{k) exp{n{k)t) (1) 

we obtain the growth rate n for every wave number k in 
our spectra, where Ai is the initial amplitude. For each fit 
it was monitored that A(k) was smaller than 40 % of its 
wavelength, the criterion for the linear regime p^ . 

In order to test the reliability of our experimental setup 
we perform 100 independent runs with one set of material 
parameters. The particular values of this set are given in 
the first row of Table I (Experiment I). We only analyze 
image series where the angular velocity of our rotating ap¬ 
paratus is larger than 6.6 rad/s. We find that 42 fast runs 
show only a slight deviation in the angular velocity: 7.4 ± 
0.1 rad/s. These 42 measurements are analyzed to obtain 
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a mean growth rate and a mean initial amplitude for each 
wave number. The results are shown in Fig. It is seen 
that the growth rates starting with small values increase 
with increasing wave number until they saturate at larger 

k. In contrast the initial amplitude decreases for increas¬ 
ing wave numbers. In the case of large wave numbers we 
do not obtain exponential fits for every experimental run. 
This is due to the fact that the amplitude is very small 
and that we approach the limit of the resolution of our 
image processing. The error bars in Fig. || indicate that 
the number of runs which can be analyzed decreases for 
larger wave numbers. 

An obvious question concerning the underlying mech¬ 
anism of the pattern formation in our system is: How do 
different material parameters effect the dispersion rela¬ 
tion? Therefore we carry out experiments with different 
material configurations which are characterized in Table 

l. Using 2 g, 4 g, and 8 g of sand (Experiment I, II, III) 
we observe a shift of all growth rates towards larger values 
with increasing mass of sand. This effect is independent 
of the wave number (Fig. ^ (a)). Using three different size 
distributions (Experiment I, IV, V) the ratio of the cell 
thickness of 4 mm to the mean particle diameter varies 
between ^ 70 and ~ 40. As figure ^ (b) shows the mean 
particle diameter does not have any significant influence 
on the dispersion relation n{k). As a common feature we 
find the same overall behavior for all material sets: The 
growth rates increase for small k and reach a plateau for 
larger values of k. 


3 Theory 


We choose a two fluid system as a model to describe the ex¬ 
perimental results. In the initial state two incompressible 
fluids of constant densities pi and p 2 and constant dynam¬ 
ical viscosities pi and p 2 are arranged in two horizontal 
strata. The index 1 (2) refers to the fluid at the bottom 
(top) of the system. The pressure is a function of the ver¬ 
tical coordinate 2 only; x and y are the coordinates in the 
plane perpendicular to z. The acceleration due to gravity 
acts in negative z direction. A linear stability analysis is 
carried out for small disturbances of this initial state. The 
instability of a planar interface z = Zs{x,y) = 0 between 
the two fluids is known as the Rayleigh-Taylor instability 


We assume that the boundaries in z direction are far 
from the interface. Small changes 5zs in the form of the 
interface cause a pressure difference which is balanced by 
the product of the surface tension Tg and the curvature 
of the interface. Considering small disturbances 5p in the 
density and Sp in the pressure the Navier-Stokes equations 
read ||^ 


pdtu = -dxSp + pAu + {dxW + dzu) dzp , (2) 
pdtv = -dySp + pAv + {dyW + dzv) dzP , (3) 
pdtw = —dz6p -b pAw + 2dzwdzP — g6p 

+Tg {dl + ( 9 ^) Szs S(z - Zs) , ( 4 ) 


where di = d/di, i = x, y, z, t. The z dependence of p 
gives rise to the third term at the right-hand side of (^ 
y) since the viscous part in the Navier-Stokes equations 
is dj [p {djVi + diVj)] for an incompressible fluid. For con¬ 
venience, we adhere to dzP though dzP is different from 
zero only at the interface. The components of the velocity 
held V are Vx = u, Vy = m and Vz = w and are consid¬ 
ered small, so that Eqs. (g]4[) contain only terms which are 
linear in the disturbances. The delta-function 5{z — Zs) en¬ 
sures that the surface tension appears at the interface Zg 
between the two fluids. The equation of continuity (mass 
conservation) for an incompressible huid is 

div V = dxU -b dyV + dzW = 0 . (5) 

Additionally the equation 

dt6p = -(vgrad)p = -wdzP (6) 

relates the temporal variations in the density huctuations 
to the density jump at the interface which moves with w in 
z direction. The Eqs. (^^ govern the linearized system. 
The disturbances are analyzed into normal modes thus 
seeking solutions which x and t dependence is proportional 
to exp(*A:a: -b nt). The wave number is denoted by k and 
n{k) is the growth rate of the corresponding mode k. If 
the fluid is confined between two rigid planes the boundary 
conditions are 

w = dzW = 0 at z = ±oo . (7) 

where we shift the planes to infinity for the sake of simplic¬ 
ity. The other boundary conditions are related to the in¬ 
terface. All three components of the velocity and the tan¬ 
gential viscous stresses must be continuous. Using the ex¬ 
ponential ansatz for the continuous velocity components, 
Eq. (H) gives the continuity of dzW, too. The continuity 
of p{d'/ + fc^) w is the condition which ensures that the 
two tangential stress components are continuous across 
the interface. Inserting the exponential ansatz into (^j^) 
and determining the solution of w in each region in such a 
way that the boundary condition ( 0 ) as well as the interfa¬ 
cial conditions are satisfied leads to the dispersion relation 




k‘^T 

L 9KPi+P2)\ 


+ 1 


(0291 + Q:i92 - k) 


— 4 ^ 010:2 H-(aii^i — Q!2J^2)[Q!29 i — Q!i92 + k{ai — 

n 

4k^ 

- 0 : 2 )] H-“ «2J^2)^(9i - k){q 2 - k) =0. (8) 

The abbreviations 01,2 = pi, 2 /{pi + P 2 ) and < 7 ^ 2 = + 
n/j/i _2 were introduced where vx ^2 = Pi, 2 /Pi ,2 is the kine¬ 
matic viscosity for each region. For the rest of the paper we 
will use a dimensionless surface tension S = T/[(pi -b P 2 ) 
{gvfy^^] if not stated otherwise. 

We now briefly discuss the general results of the dis¬ 
persion relation (@). The configuration where the lighter 
fluid is on top of the heavier one, p 2 < Pi, is always stable. 
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i.e., Re{n(fc)} < 0 for all k. This is independent of whether 
there is any surface tension (Fig. ^ (b)) or not (Fig. 0(a)). 
If the strata are in the opposite order, p 2 > pi, then the 
surface tension plavB a crucial role. In the case of no sur¬ 
face tension (Fig. g (a)) the system is unstable against 
disturbances of any wave number, i.e., Re{n(fc)} > 0 for 
all k. If there is a surface tension a critical wave number 
exists 

kc = y|;(p2 - pi) , (9) 

and the system is stable (unstable) against modes with 
wave numbers which are larger (smaller) than kc (Fig. ^ 
(b)). A moderate variation in the relation between the two 
viscosities vi and 1^2 has no strong influence on the general 
behavior of the growth rates (see Figs. 0 and ||). 

The experiments are carried out in a finite-size cell in 
contrast to our simplification of infinite length in z direc¬ 
tion. We determine now the limits to which this simplifi¬ 
cation is justified (see also [^). If the walls of the cell are 
at z = ± 1 ^ 2 : the appropriate ansatz for w is 

wi = aie+'=^ -b 020 "'=^ -b 6ie+«i^ -b for z < 0, 

( 10 ) 

W 2 = cie+'=^ -b C 2 e-^^ -b die+«i^ -b ^ 26 "'?^^ for z > 0 . 

( 11 ) 

With the boundary condition w = d^w = 0 at z = ±Lz 
and the analysis at the interface as above one ends up 
with a system of eight equations for the constants in dl- 
0) . The coefficients of the corresponding matrix are given 
in the Appendix. The vanishing determinant leads to the 
dispersion relation n{k), its numerical solution is shown in 
Fig. ^ for different lengths Lz- The comparison with the 
data for infinite Lz shows that there is no real difference 
as long as iLzl > 3 mm for the used material parameters. 
Since the thickness of the sand layer is below this margin 
we expect finite-size effects. In the case of asymmetrically 
arranged walls at z = —Li and z = +L2 the dispersion 
relation shows up to a ratio of |Li|/L 2 — 17 only a very 
weak deviation from the results for symmetrical walls. 

In the frame of our continuum approach we consider 
the water-sand mixture as a suspension in accordance with 
the classification in 1^ . The dispersion medium is water 
and the dispersed material consists of sand particles with 
a mean diameter and a density as stated in Sect. The 
material density as well as the dynamical viscosity of the 
mixture depend on the packing density (j) of the granular 
material. The packing density (j) measures the volumetric 
concentration of the particles in the mixture. The material 
density of the mixture is given by 

Pmixture — 4^Psand “b (1 Pwater • (1^) 


Since Pwater and Psand are constant, the considered small 
disturbances in the mixture density Spmixture imply that 
(j) varies according to Eq. (p^. In the following the index 
I refers to water and 2 to the mixture, respectively. Two 
empirical formulae p0|,^ 



0.63, (13) 



(t^max = 0.605,(14) 


were widely used for the dynamical viscosity of a hard 

g ’ lere suspension. The maximal packing densities in ( p^ , 
result from the fit of the proposed formulae to the 
experimental results. 

Figure ^ shows the behavior of the relative dynamical 
viscosity pr for different packing densities (j) according to 
Eqs. (|^, 1^). A third relation, ^ 1/ [1 - {4’/4>max)f^^ 
with (pmax = 0.625, was also plotted where particle sizes 
beyond the colloidal range were incorporated into the the¬ 
oretical basis of this relation All three of them give 
nearly the same value for Pr{4’) 3-f moderately dense pack¬ 
ings, 0.4 < (j) < 0.48. Above this region the relative viscos¬ 
ity diverges as cj) approaches its maximal value. But the 
strength of the divergent behavior is remarkably different 
so that pr already varies about a factor of 5 ioi (j) = 0.57. 
The fact of wide spread measured values for the relative 
viscosity of the same system is exemplarily shown in . 
The measured values of pr vary between 21 and 400 for 
a rather colloid system aX tp ~ 0.5. A similar situation is 
present for suspensions ||3^ which is why we cannot rely 
on a fixed value of pr within small error bars. Therefore we 
consider pr a variable parameter within reasonable limits 
rather than a fixed material parameter. 


4 Results and Discussion 

In experiment I the height of the sedimented sand layer 
was measured by an optical close-up with a resolution of 
37 pm per pixel. On the basis of six independent sam¬ 
ples the height was determined to 2.6 ± 0.2 mm. This re¬ 
sults in a packing densitiy of <j) = 0.48 ± 0.04. Thus the 
mixture is sufficiently characterized by Eqs. (00) where 
Pwater = 0.988 g/cm^. By means of 0) the viscosity of 
the mixture is p2 — 15/ri for (j) = 0.48. The resulting 
growth rates show disagreement for infinite Lz as well as 
for finite Lz (Fig. 0 )- Since the relative dynamical viscos¬ 
ity is the most uncertain quantity in our calculations we 
vary pr to find the best fit with the experimental data. Un¬ 
der the assumption of zero surface tension a least square fit 
results in p 2 ~ 104 pi which gives a fairly good agreement 
over the whole k range. In comparison to p2 — 104 fits 
with p2 — Pi and p2 — 85 pi show a better agreement 
for smaller and larger wave numbers k, respectively. 

For two sets of parameters, p2 — 15 and p2 — 
104/Ti, we calculate the dispersion relation for infinite as 
well as for finite Lz ■ The finite-size effects appear for small 
wave numbers and decrease with increasing wave numbers. 
The differences in n{k) for small wave numbers are more 
pronounced for larger relative viscosities (see Fig. 0 ). The 
dispersion relations for infinite and finite Lz approach each 
other at wave numbers where the initial amplitudes of the 
disturbances are in the order of 10“^ mm. For these dis¬ 
turbances a boundary at 2.6 mm distance appears to be at 
infinity. Therefore it does not matter whether we choose 
Lz = ±oo or Lz = ±2.6 mm. This is not the case for 
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small wave numbers where the initial amplitudes are in 
the order of 10“^ mm (Fig. |{b)). _ 

There is one viscosity measurement [Q which comes 
near to the value of the relative viscosity suggested by our 
fit. The measurement was carried out for crushed sand 
with particle diameters from 20 to 80 /im and gives a value 
of /r 2 ~ 110 fj,i for (j) = 0.48 (see Fig. 3 in ||3^). This value 
for the relative viscosity is notably close to our fit value. 
The values for the viscosity of the suspension estimated 
in 1^, ij ,2 ~ (130 — 190) /ii, are in a similar range as ours. 
These values are determined by the help of the singular 
wave number with the largest growth rate. A real compari¬ 
son with our fit value is not possible because the necessary 
packing density is not given for the type of experiments 
from which the viscosity values were estimated. 

Nevertheless, the large difference in between 
^ and I ^ over a wide range of the packing density (j) rep¬ 
resents an unsatisfying situation. It highlights the need for 
comprehensive and unambiguous viscosity measurements 
in highly concentrated hard sphere suspensions. It further 
shows how sensitive fir is on experimental methods |32|, 
the accurate determination of the density packing 331, 
and the type of flow involved in the measurements ||36|| . 

The aspect ratio of the cell suggests that a description 
refering to the Hele-Shaw type of the cell might be closer 
to the experimental configuration. Adapting the disper¬ 
sion relation of the Saffman-Taylor instability to the case 
of zero throughflow velocity leads to M 


n(k) = 


12(/ii + /i2) 


[kg{P2 - Pi) - Tk^] (15) 


where 6 = 4 mm denotes the width of the cell. In contrast 
to the Rayleigh-Taylor approach Eq. & contains two 
parameters, the dynamical viscosity fj ,2 of the mixture and 
the surface tension T, which have to be determined by a 
least square fit. The validity of the Hele-Shaw approach is 
limited by a cut-off condition at which the wave number 
exceeds 27r/6. 

The fit values /i 2 — 404/ri and T ~ 1.06 • 10“^ Nm“^ 
result in a ht curve which is inferior to the Rayleigh-Taylor 
ht (Fig. Additionally, both fit parameters are ques¬ 
tionable. The value ^2 — 404/ri is beyond any realistic 
one for the dynamical viscosity of the mixture at a pack¬ 
ing density of (f> = 0.48. A nonzero surface tension between 
the water-sand mixture and water is also arguable. 

The better agreement of the Rayleigh-Taylor approach 
with the experimental data is backed up by calculations fo¬ 
cused on the influence of finite cell boundaries for Rayleigh- 
Benard convection |Q. As long as the ratio between the 
relevant height of the cell and its width is smaller than 
0.8, the assumption of infinite boundaries in y direction 
is a good approximation. If this ratio is larger than 5 the 
Hele-Shaw approach is well-founded. The height of the 
sand layer and the width of the cell give a ratio of 0.65 
which supports the Rayleigh-Taylor approach. As a con¬ 
sequence the cell width has no influence on the instability 
for the used material sets. 

It has to be stressed that both the Hele-Shaw and the 
Rayleigh-Taylor approach assume a trivial y dependence 


of the flow. In the experiment, that assumption is not 
totally fulfiled: The rotation of the apparatus leads to a 
sand layer which is not perfectly flat even in the beginning 
of the flow process. This will lead to a three-dimensional 
flow which effects the appearing patterns and their wave¬ 
lengths. The strength of this effect is, however, presently 
hard to determine. 


5 Concluding Remarks 


In a closed Hele-Shaw-like cell the temporal evolution of 
a water-sand interface was investigated. For the unstable 
stratihcation, sand above water, the instability is driven 
by gravity. The images of the temporal evolution were an¬ 
alyzed by DFT. The Fourier spectra show that the initial 
disturbances of the interface grow exponentially at the be¬ 
ginning of the pattern forming process. This enables us to 
determine the growth rates by an exponential fit for ev¬ 
ery wave number in our spectra. The data show that the 
growth rate increases with increasing wave number until 
it saturates at larger values of k. This general behavior is 
not influenced by the mean particle diameter which was 
tested with three different size distributions. Experiments 
with different amounts of sand reveal a shift of the disper¬ 
sion relation n{k) as a whole. By increasing the mass of 
sand, n{k) is shifted towards larger values of n. 

To describe the general behavior we choose a two-fluid 
system as a model. Carrying out a linear stability analysis 
for the interface between the two fluids we calculate the 
growth rates from the dispersion relation for a finite-size 
cell. The theoretical results agree with the essence in the 
experimental findings when assuming a relative viscosity 
of the water-sand mixture which is close to one measured 
with crushed sand and water |^ . 

Considering our simplihcations and the uncertainty in 
one relevant material parameter, the continuum approach 
gives a reasonable agreement with the experimental re¬ 
sults. An interpretation for such an agreement was pro¬ 
posed in . Above a critical value of the packing density 
the mixture exhibits a non-zero yield stress well known for 
dry granular material. Below this critical value the yield 
stress vanishes and the mixture behaves more like a fluid 
with effective properties. Applying this interpretation, the 
measured packing density of (f> = 0.48 is below this critical 
value. Thus, our results confirm a certain analogy between 
concentrated suspensions and fluids which was found also 
in numerical simulations |24 . However, further investiga¬ 
tions need to be done in order to clarify open questions. 

The available experimental data give no hint on whether 
the dispersion relation n{k) will have a second zero at 
kc after the observed plateau. A critical wave number kc 
means a non-zero interfacial tension according to (^). A 
surface tension acts in a way to minimize the surface of 
the fluid. Therefore it suppresses the formation of waves 
with large wave numbers because their creation entail ad¬ 
ditional surface. The suppression leads to a saturation (re¬ 
duction) in the growth rates of disturbances with large 
(very large) wave numbers. In the context of miscible flu¬ 
ids with slow diffusion the concept of an effective dynam- 
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ical surface tension was recently successfully applied . 
The reason for such a surface tension lies in the attraction 
between moving particles in a fluid for non-zero Reynolds 
numbers ||40|| . The attraction originates from the dynamics 
in the viscous fluid and results in a favoured distance be¬ 
tween the moving particles. Consequently, is costs energy 
to separate the particles beyond this favoured distance, 
i.e., to dilute the suspension. The necessary energy corre¬ 
sponds to a surface tension which is called effective dynam¬ 
ical surface tension to emphasize its dynamical origin. The 
effect of such an effective dynamical surface tension is seen 
in experiments by the presence of sharp interfaces in rising 
bubbles, falling drops M and of growing deposits 39 . A 


concentration gradient concerning sand is present in our 
system between the water and the water-sand layer. By 
reducing the density difference between the layers, mea¬ 
surements for k values beyond the plateau should be pos¬ 
sible. This may lead towards the determination of kc and 
the interfacial tension, respectively. 

Since a linear stability analysis is restricted to terms 
linear in the disturbances we cannot take into account 
fluctuations in the relative viscosity of the mixture. These 
would lead to terms of higher order because all terms 
which contain the viscosity are already linear in the veloc¬ 
ity disturbances (see |[^. On the other hand, considering 
fluctuations of Pmixture(4>) ^ huctuations of the packing 
density 0, implies variations of This inconsistency 

owing to the restrictions of our linear theory can only be 
resolved in a nonlinear analysis. 


We are grateful to Lluis Carrillo, Stefan Neser, and Stefan 
Schwarzer for inspiring discussions. The experiments were sup¬ 
ported by DFG through Re 588/11-1. 
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Neon tube 




Fig. 2. Sand-water interfaces at certain time steps: (a) 20 ms, 
(b) 80 ms, (c) 140 ms, (d) 200 ms, and (e) 260 ms. The size 
distribution of the sand particles is given by 71-80 ^m and the 
sand mass in the cell by 2 g, respectively. The presented frames 
show the middle part of the cell and have a horizontal length 
of 68 mm. 


Fig. 1. Experimental setup. 
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Fig. 3. Temporal evolution of the water-sand interfaces. The 
patterns are detected every 20 ms and shown with a constant 
vertical offset of 1 mm. The experimental conditions are the 
same as in Fig. ^ 



Fig. 4. Amplitude A of a DFT-analysis for a typical wave 
number (here k = 7 cm“^) in dmendence on time t. An ex¬ 
ponential fit is obtained by Eq. (|]l). The values belong to the 
interfaces presented in Fig. ^ 



Fig. 5. Growth rate (a) and initial amplitude (b) versus the 
wave number k for experiment I (see Table 
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Fig. 6. Growth rate versus the wave number for different ma¬ 
terial parameters. In (a) 2 g (□), 4 g (•), and 8 g (a) of sand 
with the same size distribution were used (Experiment I, II, 
III). In (b) 2 g of sand with size distributions of 56 — 63 ^m 
(□), 71 — 80 pm (a), and 90 — 100 pm (•) were used (Experi¬ 
ment IV, I, V). 


Fig. 7. For zero surface tension (a) and surface tension 5=1 
(b) the growth rate Re{n} of the modes is plotted versus their 
wave number k for the case that the top fluid is lighter than 
the bottom one, p 2 = 0.5pi. The surface tension as well as the 
relation of the viscosities does not change the overall behavior 
too much as the three examples show: v\ = V 2 = 10“®m^s“^ 
(solid line), vi = 10“®m^s“^ and V 2 = 0.75zzi (long-dashed 
line), and V 2 = 10“® m^ s”*^ and = 0.75iZ2 (dot-dashed line). 
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Fig. 8. The growth rate Re{n} of the modes is plotted versus 
their wave number k for the case that the top fluid is heavier 
than the bottom one, p 2 ~ 2pi. The surface tension S is zero 
in (a) and 1 in (b). The surface tension causes a drastic change 
in the behavior whereas the ratio of the viscosities does not 
change the overall behavior too much as the three examples 
show: = V 2 = 10“®m^s“'^ (solid line), v\ = 10“®m^s“'^ 

and V 2 ~ 0.75izi (long-dashed line), and U 2 — 10“®m^s“^ and 
= 0.75 i^ 2 (dot-dashed line). 
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Fig. 9. For zero surface tension (a) and surface tension S = 
12.8 (b) the growth rate Re{n} is plotted versus the wave num¬ 
ber k for Lz — ±3 mm (long-dashed line) and Lz — ±lmm 
(dashed line). For comparison the solid line shows the graph 
for inhnite Lz- It reveals that finite-size effects play a signifi¬ 
cant role if fc < 27r/|L2j. The material parameters are p 2 = 2pi 
and V 2 = 33izi with pi = lgcm“® and vi = 10“® m^ s“^. 
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Fig. 10. Relative dynamical viscosity versus the packing 
density cj). The results of the two empirical formulae (|^, p^ ) are 
plotted as dashed and long-dashed lines, respectively. The solid 
line indicates Hr ^ 1/[I — with ij)max = 0.625 

I p^ . Hr is nearly the same for all three approaches provided the 
packing is not too dense, (j) < 0.48. Above this range Hr starts 
to diverge as (j) reaches (jimax where the divergent behavior 
differs signihcantly between the various approximations. 



Fig. 11. Growth rates n of the disturbances against their wave 
numbers k for a packing density cj) = 0.48. The results n{k) for 
infinite Lz according to are shown by the dashed line [Hr = 
15) and the long-dashed line {Hr = 104), respectively. The 
solid lines are calculated with \Lz\ = ±2.6mm, the measured 
height of the sand layer. The relative viscosity for each drawn 
curve is stated at the right end of it. The finite-size effects 
appear for small wave numbers and increase with increa sing 
relative viscosity. A relative viscosity of Hr = 15 based on (p^ 
gives growth rates which are far away from the experimental 
results (■). The best fit over the whole k range gives Hr ~ 104. 
Fits with slightly smaller or larger Hr values deliver better 
agreements either with larger or smaller wave numbers. 



Fig. 12. Growth rates n of the disturbances against their wave 
numbers k for a packing density cj) = 0.48. Two different fits 
are compared with the experimental results (■). The solid line 
shows the graph based on the Rayleigh-Taylor approach with 
finite z boundaries at \Lz \ = ±2.6 mm and Hr = 104. The long- 
dashed line presents the results of the Hele-Shaw approach with 
Hr = 404 and T ~ 1.06 • 10“® Nm'^ 








